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Abstract 

We survey the equations of continuous-time quantum walks on simple one-dimensional lattices, 
which include the finite and infinite lines and the finite cycle, and compare them with the classical 
continuous-time Markov chains. The focus of our expository article is on analyzing these processes 
using the Laplace transform on the stochastic recurrences. The resulting time evolution equations, clas- 
sical versus quantum, are strikingly similar in form, although dissimilar in behavior. We also provide 
comparisons with analyses performed using spectral methods. 

1 Introduction 

The theory of Markov chains on countable structures is an important area in mathematics and physics. A 
quantum analogue of continuous-time Markov chains on the infinite line is well-studied in physics (for 
example, it can be found in 1 12 1, Chapters 13 and 16). More recently, it was studied by Aharonov et al. |2| 
and by Farhi and Gutmann 1 1 1 1. The latter work placed the problem in the context of quantum algorithms for 
search problems on graphs. Here the symmetric stochastic matrix of the graph is viewed as a Hamiltonian 
of the quantum process. Using Schrodinger's equation with this Hamiltonian, we obtain a quantum walk on 
the underlying graph, instead of a classical random walk. 

Recent works on continuous-time quantum walks on finite graphs include the analyses of mixing and 
hitting times on the n-cube II16II14L of mixing times on circulant graphs and Cayley graphs of the symmetric 
group l6l fT3ll . and of hitting times on path-like graphs |8l|9l. Most of these are structural results based on 
spectral analysis of the underlying graphs, such as the n-cube, circulant and Cayley graphs, and (weighted) 
paths. For example, Moore and Russell 1 16 1 proved that the mixing time of a quantum walk on the n-cube 
is asymptotically faster than a classical random walk; Kempe llT4l proved that the hitting time for vertices 
on opposite ends of the n-cube is exponentially faster than in a classical random walk. Ahmadi et al. O 
and Gerhardt and Watrous fT3l proved that circulants and the Cayley graph of the symmetric group lack the 
uniform mixing property found in classical random walks. 

A recent work of Childs et al. [9| gave intriguing evidence that continuous-time quantum walk is 
a powerful method for designing new quantum algorithms. They analyzed diffusion processes on one- 
dimensional structures (finite path and infinite line) using spectral methods. Another work by Childs and 
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Goldstone fTOl explored the application of continuous-time quantum walks to perform Grover search on 
spatial lattices. 

There is an alternate theory of discrete quantum walks on graphs, which we will not discuss here. This 
alternate model was studied in Aharonov et al. and Ambainis et al. 0, but had appeared earlier in work 
by Meyer fTSl . The work by Ambainis et al. 13 had also focused on one-dimensional lattices. Recently, 
Ambainis |T| developed an optimal (discrete) quantum walk algorithm for the fundamental problem of 
Element Distinctness. This offers another idea for developing quantum algorithms. 

We survey and (re)derive equations for the continuous-time classical and quantum walks on one- 
dimensional lattices using the Laplace transform that works directly with the recurrences. The Laplace 
transform is a well-known tool in stochastic processes (see Q) and it might offer a useful alternative to the 
Fourier transform in certain settings. 

1.1 Stochastic walks on graphs 

Let G = {V, E) be a simple (no self-loops), countable, undirected graph with adjacency matrix A. Let D 
be a diagonal matrix whose j-th entry is the degree of the j-th vertex of G. The Laplacian of G is defined 
H = A — D. Suppose that P{t) is a time-dependent probability distribution of a stochastic (particle) 
process on G. The classical evolution of the continuous-time walk is given by the Kolmogorov equation 

P'{t) = HP{t). (1) 

The solution this equation, modulo some conditions, is P{t) = e*^P(0), which can be solved by diagonal- 
izing the symmetric matrix H. This spectral approach requires full knowledge of the spectrum of H. 

A quantum analogue of this classical walk uses the Schrodinger equation in place of the Kolmogorov 
equation. Let ip : V{G) ^ C be the time-independent amplitude of the quantum process on G. Then, the 
wave evolution is 

i^j^m = Hm- (2) 

Assuming /i = 1 for simplicity, the solution of this equation is ip{t) = e~**^V'(0), which, again, is solvable 
via spectral techniques. The classical behavior of this quantum process is given by the probability distribu- 
tion P{t) whose j-th entry is Pj{t) = \ipj{t)\'^, where t/^jit) = {j\'ijj{t)). The average probability of vertex 
j is defined as P{j) = Huit^oo ^ Jq Pj{t)dt (see |4|). 
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Figure 1: Examples of some one-dimensional lattices. From left to right: Z, P4, Z4. 



The table in Figure |2l shows the known equations for continuous-time stochastic walks on the infinite 
(integer) line Z = {. . . , —2, —1, 0, 1,2,.. .}, the finite cycle Z^r = {0, . . . , — 1} on iV vertices, and the 
finite path P^r = {0, . . . , N} on + 1 vertices, in terms of the two kinds of Bessel functions /(•) and J(-). 
We assume here that the particle is initially at 0. The plots in Figures |3l and |4] show the dissimilar behavior 
of the classical versus quantum walks. 
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Graph 


Classical walk 

Pj (t) = probability on vertex j at time t 


Quantum walk 

ipj (t) = amplitude on vertex j at time t 
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Figure 2: The equations of the continuous-time classical versus quantum walks on the infinite line, finite 
cycle, and the finite line, assuming the particle starts at position 0. 

1.2 Laplace transform 

The Laplace transform of a time-dependent function Pit), denoted P(s) = is defined as 

/■oo 

£{P(t)} = / e-"*P(t) dt. (3) 

^0 

The only basic properties of the Laplace transform which we will need are (see 0): 

• Linearity: L{aP{t) + = aP(s) + hQ{s) 

• Derivative: L{P' (f)] = sP{s) - P{0) 

• Shifting: £{e"*P(t)} = P(s - a) 

The relevant Inverse Laplace transform involving the Bessel functions are (for u > —1): 

p^^^ ^ {s - Vs^ - a ^ f_ ^ ^^^^ ^ a^I^iat) (Eqn. 29.3.59 in H) (4) 

vs^ — 

p^^^ ^ Ws^ + a? - s Y_ ^ p^^^ ^ a" Mat) (Eqn. 29.3.56 in 0) (5) 

2 The infinite line 

Classical process. The Kolmogorov equation for the infinite line is 

= - P,{t) + (6) 

with initial value Pj{0) = Sqj. The Laplace transform of Q is 

Pj+i{s) - 2(8 + l)Pj{s) + = -Pj{0). (7) 
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The solution of — 2(s + \)q + 1 is g-t = (s + 1) ± ^/{s + \Y^^. A natural guess of the solution is 

p.ls)-{ M ifi<0 (8) 

When j = 0, we get A = {1 + s — q-)~^. Thus, for j £ Z, 

p . . _ g^'' _ ((. + i)-V(. + 1)^-1)1^-1 

Using the Inverse Laplace transform @, after shifting 5 = s + 1, we get 

P,(t) = e-%|(t). (10) 
This is a probability function, since e*/2(^+^/^) = Yl'h=-oo ^^^k{^)> if 2; / (see Eqn. 9.6.33 in l|3J). 

Quantum process. The Schrodinger equation for the infinite line is 

iV-W = \iljj-i{t) + \^j+i{t). (11) 

The Laplace transform of (fTTt is 

- 2i{sijj{s) - V'j(O)) + ^j^i{s) = (12) 
The solutions of q^ — 2isq + 1 are g± = i{s ± V s'^ + 1), where q+q- = 1. A guess for the solution is 

V- ifj>0 ^ ' 

When j = 0, we get A = {s + iq )^^. Thus, 

The Inverse Laplace transform (ISjl yields, for j S Z, 

V',-(t) = H)i^V|,-|(t), (15) 

This is a probabiUty function, since 1 = J^{z) + 2^^^^ J|(z) (see Eqn. 9.1.76 in H). 

Spectral analysis. Let i7 be a Hamiltonian defined as {j\H\k) = ^ if j = k ± 1, and otherwise. For 
each p £ [— vr, vr], define \p) so that 

{j\p) = -^e'P^. (16) 
v27r 



The eigenvalue equation H\p) = Xp\p) has the solution Xp = cos(p). Thus, the amplitude of position j 
when the particle starts at position is 

{j\e-'^'\0) = — r e^fJe-^*™^(f) dp = {-iyjj{t) (see Eqn. 9.1.21 in 0) (17) 

271" J-n 

Childs et al. O gave a more generalized analysis along these lines. 



4 



3 The finite cycle 



Classical process. If A is the adjacency matrix of the finite cycle, let H = — I he its Laplacian matrix. 
The Kolmogorov equation for the finite cycle is 



Pj(t) = iP,_i(t)-P,(t) + iP,+i(t). 
Applying the Laplace transform to (fisl . we get 



(18) 



(19) 



For convenience, define the extra condition = P/v^i(s) + 2, so that Pj_|_i(s) — 2(s + l)Pj{s) + 

Pj-i{s) = holds for j G Zjv- The cycle condition is Pn{s) = Po{s). We guess the solution to be 



Pj{s)=Aq'^+Bq>_, 



(20) 



where q± is the solution to — 2(s + \)x + 1 = 0, i.e., q± = {s + 1) ^ + 1)^ — 1> with g+g- = 1. 
Using the cycle condition, we get 



Aq^ + Bq'l = A + B 



B = Aq^^. 



> A{q^-l)=B{l-q'^) = 
Using the extra condition and SITh . we get A = 2{{q^ — q^){q^ — 1))^^. Thus, for j G Zat, 



(21) 



Pj{s) = Aq\_ + Bqi_ = A{ql + q 



2 (d + q-~'] 



fc=0 



((s + 1) - V(s + 1)2 - If^+i {{s + l)- v^(s + l)2-l)(fc+i)^-i' 



+ 1)2-1 



+ 



The Inverse Laplace transform Q, after shifting, yields, for j G 



Pj ) = E t^'^^+J- + hk+l)N~j (*)] = 5^ (t) . 

fe=0 Q=±j(modAf) 



(22) 



Quantum process. Since the finite cycle is a regular graph, instead of the Laplacian, we use the adjacency 
matrix directly. In a continuous-time quantum walk, this simply introduces an irrelevant phase factor in the 
final expression. The Schrodinger equation, in this case, is 

^■(t) = \^j-i{t) + l^j+i{t). (23) 

The Laplace transform of d23b is 

^j+i{s) - 2i{si^j{s) - V'j(O)) + V^j-i(s) = (24) 
The cycle boundary condition is iIjn{s) = V'o('S)- For convenience, define 

i^^i{s) = 'ipN~i{s) = 2i. (25) 
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The solutions of — 2isq -\- 1 are q± = i{s ± + 1), with q+q^ = 1. A solution guess, for j G Zjv, is 

iPj{s)=Aq{ + Bqt. (26) 

The cycle boundary condition yields B = Aq'^. By (l25l . we get A = 2i(((7_|_ — q-){q^ — 1))"^- Thus, for 

j G Zn, 



{s) = Aq\ + Bql = A{ql + q" 



2i iqi_ + q 



2i 



iq+-q-) (i-q^) {q+-q-) 



k=0 



E 

fe=0 



+ 1 + 1 

The Inverse Laplace transform ^ gives, for j G Zn, 

oo 



fc=0 



E i-iTJait). (27) 

Q=±j(modiV) 



Spectral analysis. The normalized adjacency matrix of Z^r is the circulant matrix 
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(28) 



It is well-known that all A'^ x A'^ circulant matrices are unitarily diagonalized by the Fourier matrix F = 
-^V{ijJn), where ujn = e^'^*/^ and V{ijJn) is the Vandermonde matrix defined as V{ijJN)[j,k] = w^, 

for j,k G {0, 1, ... — 1}. The eigenvalues of H are \j = ^(wj^ + oj^^^ = cos{2TTj /N), for 
j = 0, 1, . . . , AT — 1. Thus, the wave amplitude at vertex j at time t is 



N-l 



AT 



~itcos{2iTk/N) jk 



N ■ 



(29) 



fc=0 



From earUer analysis, we get the following Bessel equation 



N ^ 

k=0 



^~it cos(2-Kk/N) ^2-Kijk/N 



o=itj (mod TV) 



(30) 



It is an open question if there exists a time t G such that for all j G Z^ we have = l/-^» i-C-. 

uniformity is achieved at some time t. For = 2, 3, 4, it is known that instantaneous exact uniform mixing 
is achieved (see If6l l6l'). 
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4 The finite path 

Classical process. Let A be the normalized adjacency matrix of the finite path, where A is a stochastic 
matrix with the probability transitions are proportional to the degrees of the vertices. Let H = A — Ibeits 
Laplacian. Then, the Kolmogorov equation, in this case, is 

p;it) = - Pj{t) + ip,+i(t), (31) 

for < j < A^, with initial condition Pj{0) = 5j,o and boundary conditions 

P^(t) = P,{t) - Po{t), Pl^it) = PN-i{t) - PN{t). (32) 
The Laplace transform of d3l1) is 

P,+i(s)-2(s + l)P,(s)+P,„i(s) = 0, 0<j<N, (33) 

and two boundary equations (1 + s)Po{s) — 1 = Pi{s), and (1 + s)Pn{s) = P/v-i(s)- A guess of the 
solution is 

Pj (s) = Aq{_ + Bqi , < j < A^, (34) 



where q± = (s + 1) ± y'Js + T)^^^. The boundary equations give B — A = 2/{q^ — q^) and A = Bq^ . 
Combining these last two equations, we get 

2 ql^ 

^ - (g+-g_)(l-e^)- ^^^^ 

Thus, for j = 0,1,..., AT, 

P,{s) = Aqi^+Bc^_=AW^ + qf-^^- ^ i^'- + I-""'') 



oo 



2 ^^^2Arfc+i ^ ^2N{k+l)-j^ 



The Inverse Laplace transform Q, after shifting, yields, for j = 0, 1, . . . , A'^, 

oo 

Pj it) = Y.^^' [^2JVfc+, {t) + hN{k+l)-j (*)] = Yl • (36) 

fc=0 a=±j{moi2N) 

Quantum process. The Schrodinger equation for the finite path is 

iV'-W = iVi-iW + ^V'i+iW, (37) 
for < j < N, with initial condition ipj{0) = 6oj and boundary conditions 

iV'o(t) =V'i(t), i^lj'Nit) = i^N-i{t). (38) 

The Laplace transform of d37b is 

i>j+i{s)-2isil^j{s)+i>j^i{s) = 0, 0<j<N, (39) 
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and two boundary equations istl^o{s) — i = i^iis), and isijjNis) = t/^N-iis)- The solutions of — 2isq + 1 
are q± = i{s ± \/s^~+T). A guess of the solution is 

ijj (s) = Aq{ + Bqi_ , 0<j<N. (40) 

From the boundary equations, we get S — ^ = 2i/{q^ — q^) and B = Aq"^ . Thus, 



A 



2i 



[I - qi' 



2N-) 



Fori = 0,1,... ,iV, 



i;j{s) = Aq^ + Bql = A{qi + qX 



2N-j. 



2i {qi_ + 



2i 



°° \ If A\/ / , 1 „\\2Nk+i II - 

E 



^ ^^2Nk+j ^ Jik+l)N-j 



k=0 



The Inverse Laplace transform ^ yields, for j = 0, 1, . . . , A^, 

oo 

V-.C*) = E [i-i?'''^^J2Nk+j{t) + {-i?''^'^'^-^J2N(k+l)-j{t) 
k=0 



(41) 



{-iTUt). (42) 

a=±j(mod 27V) 



Spectral analysis. The spectrum of a path on n vertices is given by Spitzer flTl . For j G {0, 1, . . . , A^}, the 

eigenvalue \j and its eigenvector Vj are given by 



,'(i + l)7r\ 
Xi = cos ( —zz — — I , VjW 



N + 2 J ' ' \ N + 2 
The probability of measuring vertex at time t is given by 



{N + 2) 



j,k 



sm 



2 ^ (i + l)7r 



sm 



iV + 2 



+ l)7r 



A^ + 2 



-it{Xj-Xk) 



Since all eigenvalues are distinct, the average probability of measuring the starting vertex is 



P(0) 



{N+2rj^ 



sm — sm I — I iim 



(W + 2)2 



E 



sm 



N + 2 
A^ + 2 



AT + 2 y T^oo r Jo 



T 



-it(Aj-Afe) 



Equating this with (I42t . we obtain a Bessel-like equation: 



lim — 



a=0{mod 2Af) 



6t 



{N + 2f 



N 

E 

A:=0 



sm 



4 + l)7r 
iV + 2 



(43) 



(44) 



(45) 
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5 Conclusions 



This expository survey reviews equations for the continuous-time quantum walks on one-dimensional lat- 
tices. The focus was on analysis based on the Laplace transform which works directly with the stochastic re- 
currences. It would be interesting to extend this analysis to higher-dimensional or to regular graph-theoretic 
settings. Another interesting direction is to consider lattices with defects and weighted graphs |9|. 
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Figure 3: Stochastic walks on the infinite line Z: (a) plot of Po{t) in the continuous-time random walks for 
t G [0, 50]. (b) plot of |V'o(OP ^ continuous-time quantum walk for t G [0, 50]. Both processes exhibit 
exponential decay, but with the quantum walk showing an oscillatory behavior. 
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Figure 4: Stochastic walks on the finite cycle Z5, each approximated using 21 terms: (a) plot of Po{t) in the 
continuous-time random walks for t G [0, 50]. (b) plot of |^/;o(OP i'^ the continuous-time quantum walk for 
t G [0, 500]. The classical walk settles quickly to 1/5, while the quantum walk exhibit a short-term chaotic 
behavior and a long-term oscillatory behavior below 0.1. 
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